Resource Handbook $.50 


DIGITAL EQUIPMENT CORPORATION 


HUNTINGTON Il Simulation Program — CHARGE 





os are for kids dlilg|i|tlall —EduSystems—expandable, economical 


DIGITAL EQUIPMENT CORPORATION 
EDUCATIONAL PUBLICATIONS 


A partial list of the publications in the continuing series of curric- 
ulum material published by Digital for use with EduSystems and RSTS 
are listed below. Please inquire directly for prices on classroom 
quantities. 


Additional publications may be obtained from: 
Software Distribution Center 


Digital Equipment Corporation 
Maynard, Massachusetts 01754 





Populution: Self teaching BASIC Workbook $2.00 

BASIC Matrix Operations, Project Solo $1.00 

Problems for Computer Mathematics Shs 20 

Advanced Problems for Computer Mathematics 2.00 

BASIC Application Programs | - 
Mathematics I $1.00 
Mathematics II 1.00 
Science 1.00 
Business and Social Studies 1.00 
Plotting 1.00 


Huntington I Application Programs 
MATHEMATICS $2.00 
TEACHER ASSISTANCE LOG 


Huntington I Simulation Programs 





BIOLOGY $1.00 
CHEMISTRY 2.00 
EARTH SCIENCE 1.00 
PHYSICS 2.00 
SOCIAL STUDIES | 1.00 


Huntington II Simulation Modules 


Student Workbook $0.30 
Teacher's Guide -o0 
Resource Handbook -50 


Program Paper Tape : 1.00 


Sy 








(| CHARGE : 
YD )) 2 
| uf 
oO 
SNK Y 
YY 
YAY 
ZGI 
Y, Sal 
Le Pe 
J, YZ, NK 
Z NA 
Ay 
BUNS 
A Ny 
| ey NT” Wh 
4 ee bya g 8 
| Te NY ean 
TRG EER <S 2SRRA Ss, 
KF 40% tie 
gay wr # ) 
esa fea oe 
> RESOURCE HANDBOOK 
. 
HUNTINGTON TWO COMPUTER PROJECT 
Copyright C 1971, State University of New York 
1 September 1971 
The work of the Huntington Two Computer Project is partially 
supported by the National Science Foundation, Grant GW-5883 
. 4 


Exclusive publishing rights granted to Digital Equipment Corporation. 
REPRODUCTION NOT PERMITTED 


CHARGE 


APPENDICES 


Millikan's Oil Drop Experiment: 


The CHARGE Model .. 
Program Listing ose 
Sample Runs .... . 


Detailed Calculations 
Sg Se ee 


L&VGr 2 vss Ss 


History . 





-12 


~14 





¢ 


@r- 


APPENDIX A. Millikan's Oil Drop Experiment: History 


The results of J.J. Thomson's electron experiments (1897) and the Ruther- 
ford scattering experiments (1909) led to a more general acceptance of 
the "atomic structure of matter." During the same period many attempts 
were being made to determine whether electricity had the same proper- 
ties; that is, did there exist a basic unit of electricity? One of 

the earliest attempts along this line of inquiry was a series of exper- 
iments performed at the Cavendish Laboratory in England by J.J. Thomson. 


Just prior to Thomson's investigations, C.T.R. Wilson (also of the 
Cavendish Laboratory) discovered that air saturated with water vapor will 
begin condensation on ions introduced into the system. Starting with 

a saturated mixture,Wilson would ionize air particles by beaming x-rays 
into the solution. Shortly thereafter, he would expand the volume to 
create a supersaturated solution which would then condense on the ions 
created. The complete unit including the volume-expansion mechanism is 
known as the Wilson "cloud chamber." : 


Using a cloud chamber, Thomson noted the movements of the droplets 
formed as they drifted downward under the influence of gravity. He 
was able to determine the total quantity of water involved in the 
condensation process by measuring the amounts of water needed to sat- 
urate the volume before and after expansion. In addition, Thomson was 
also able to determine the average size of the droplets by applying 
a method devised by G.G. Stokes involving the drift of spheres (the 
droplets) through a viscous medium (the air). Thus, the total number 
of droplets produced could be determined by dividing the total quantity 
of water condensed by the mean volume of the droplets formed. By then 
introducing an electric field in the region, Thomson could separate and 
collect the droplets with the different charges to determine the total 
negative or positive charge carried by the drops. These experiments 
permitted Thomson to determine the average number of charges per drop. 


Similar experiments were performed by H.A. Wilson who used a pair of horiz- 
ontal parallel plates in the cloud chamber to produce a uniform electric 
field. The field was directed so that it opposed the gravitational 

force on the negative ions. By varying the voltage applied to the plates, 
H.A. Wilson was nearly able to balance the drop between the electric 

and gravitational forces. Determining the average size and mass of 

the drops as before, he was able to estimate the charge carried by the 
drops. 


In 1911, R.A. Millikan improved Wilson's method by. including a precedure 
for measuring the movements of a single droplet over extended periods of 
time. Also, noting that the water droplets evaporated rapidly, he 


designed his experiments to work with oil instead of water. The device 
used by Millikan is similar to that used in the modern version of the 
experiment. | 


Oil from an atomizer was sprayed above the upper plate of the apparatus 
and some of the oil drops fell through the hole to the region between 
the plates. The size of the drops could be determined by noting their 
rate of fall through air, and the mass of each drop determined from the 
size of the drop and the density of the oil. Once produced, the drops 
could be ionized by x-rays and prevented fron falling out of view of 

the observing windows by applying the proper voltage to the metal $ 
plates. Actually, when the electric fieldwas applied, some droplets 
would tend to drift upward while others increased their rate of fall. 
The speed at which the drops moved depended on the charge they carried 
and the magnitude of the field. By increasing the voltage applied to 
the plates, Millikan could sweep the faster moving particles out of 

view leaving the slower moving particles behind for further observation. 
Of the drops remaining, Millikan was able to single out one drop and 
watch it rise and fall under the influence of the various forces. 





In the actual experiments, Millikan watched single drops for hours 

at a time as they rose and fell between two cross-hairs of known 
separation marked off on his observation telescope. By noting the 

time required to move across the fixed distance over and over again, 
Millikan established "means" for the terminal velocity in each direction. Ww 
Repeating the procedure with many different single oil drops, Millikan 
obtained values for the charges that fell in well-defined groups; 

all the charges were integral multiples of a basic quantity. The 
interpretation of such results was that electric charge is not distributed 
uniformly, but comes in "clumps" or atoms of electricity" as first 
postulated by Michael Faraday fifty years earlier. | 


Once the supporting evidence for a basic electric charge had been es-_ 
tablished, determination of the number to be assigned to the basic charge 
could also be found from Millikan's experiment. The terminal velocity under 
free fall is given in terms of Stokes' constant K and the resultant 

force acting on the drop. The resultant force F is known to be: 


¥ 

Lae mg [1] 4 
where m is the mass of the oil drop. Since the density of the oil (d) © 
is known, the force F may also be written as: . : 

F = 4/3 IWRdg [2] 
With Stokes' constant K = 67ThR, 

3 
c / 677 hR | [3] 








Solving for R yields: 


ee ee [4] 


Once the radius of a given drop is known, Stokes' constant can be cal- 
culated and, subsequently, the "basic electric charge." 


One of the necessary measurements in Millikan's original experiment was 
the motion of the oil drop as it fell freely through the air. This meas- 
urement permitted the size of the drop and, consequently, the mass of 

the drop to be calculated. This additional complication can now be 
eliminated by using tiny latex spheres of known mass and size. (These 
latex spheres are ordinarily used to calibrate the magnifying power of 
electron microscopes). The spheres that are commonly available have 

a diameter of one micron (1 x 10°° meters) plus or minus about 5%, 

and a density that is very close to 1000 kg/m~. The mass of such a 
sphere is about 5 x 10-16 kilograms. The stability, uniformity and known 
radius of these spheres make possible a laboratory experiment that is 
simpler conceptually and easier to do than Millikan's original version. 
It is this experiment that is simulated in the CHARGE program. 





APPENDIX B. The CHARGE Model 


I. List of Variables and Constants 


A. Variables BASIC Notation 

1. v = velocity of the drop in meters/sec -- a 
= 
2. V = voltage applied in volts (between | V : 
-1000 and +1000) 
© 
3. c = charge on a drop in coulombs -- 
4. N = number of unbalanced electron N(1) 
charges on a drop for drop I 

B. Constants 

io. = 1.5 = 10> newton-sec/m* viscosity of air H 
at 18°C 
2 
2. g = 9.8 meter/sec acceleration of G 
| gravity 
3. q=1.6 x 10719 coulomb electron charge Q = 
4. d= 1.0 x 10° kg/m? density of latex D 
5; Rh eS es 107/ meters + 1% radius of drop R(I) 
6.D=2x ee meters plate separation _ DI 
Ii. Equations 

A. Drop velocity as a function of voltage 
1. v = F/K on 
2. oF = 7%, F, = cE r= " 
3. K = 67FhR 
4. c= Nq 
5. E = V/D 


6 a = 4/3 - 9 








| B. Charge calculation at voltage at which v = 0 
| * cE = mg 


c = mg/E = mgD/V . 


III. Number of Charges per Drop 


The number N of charges per drop is a normally distributed random 
variable with mean and standard deviation both equal to 3. 
It is determined by using the INT function on the following expression: 


3+ 3[cos(27r,)] . -21n,r5 (1) 


where r and ry are generated by the RND function. 


The function 


x = [cos(2 Wr, )] 4 -21n,r, 7 (2) 


yields very close to a normal distribution with mean equal to zero 

and standard deviation equal to 1; with the transformation 

y = ox +y , where O is the standard deviation and U the mean, 
a 7 we obtain (1). The limits on N in this program are 0< N<8. 





IV. Statistical Variation in Drop Radius 





Available latex spheres have a spread in radius from sphere to 
sphere of about + 5%. Since the cube of the radius enters into 
the equation for the mass of the sphere, the masses can vary by 
+154. In the calculation of the charge on a stopped drop, the 
average mass of a drop is used, since the individual mass is not 
known. The 154 uncertainty in mass leads to a 15% uncertainty in 
calculated charge, which masks the discreteness of charge when 
a sphere contains many charges. (For instance, for a sphere. 
S containing 8 charges, the calculated charge can be anywhere | 
&. | between 6,8 (8 - 15%) and 9,2 ( 8+ 15%). To avoid this problem 
| in the CHARGE program, the radii are assumed to have a variation 
of only + 1%. In the program the radii are uniformly distributed 
@ about the average radius R, by calculating the radius of the Ith 
drop as: : 


- 
~ 


R(I) = R(.99 + .02 RND(I)) 


where R, = 5 x 10-/ meters. The differences in radii lead to a 
spread ?n drop mass, terminal velocity, and calculated charge. 











pe ; 


199 REM NCI)=NUMBER OF ELECTRONIC CHARGE UNITS FOR PARTICLE I 

195 REM COPYRIGHT 1971 - STATE UNIVERSITY OF NEW YORK 

115 REM (NORMALLY DISTRIBUTED WITH MEAN 3 AND DEVIATION 3) 

120 REM MCI)=MASS OF PARTICLE I» BASED ON RCI) 

125 REM RCI)=RADIUS OF PARTICLE I | 

130 REM CAVERAGE RADIUS R»s 1 PERCENT UNIFORMLY DISTRIBUTED ERROR) 

135 REM MODEL DEVELOPED BY Ae CAGGIANO AND De SCARL | . 
14G REM PROGRAMMED RY C. LOSIKs JULY 1971 : 
145 REM LATEST REVISION: 8-25-72 | 

156 REM Q=RND(-1) 

155 RANDOMIZE 

169 DIM N(C4)5M(4)sRC4) 

165 PRINT " ","MILLIKAN OIL DROP EXPERIMENT" 

170 PRINT 

175 REM V:VOLTAGEs INITIALLY ZERO 

18H LET V=9 | 

185 REM H:VISCOSITY COEFFICIENT BETWEEN LATEX SPHERE AND AIR 
190 LET H=-990918 

195 REM G:ACCELERATION DUE TO GRAVITY 

200 LET G=9.8 | 

295 REM Q:COULOMBS PER ELECTRONIC CHARGE UNIT 

210 LET Q=1+6E-19 | 

215 REM D:DENSITY OF LATEX SPHERE 

2°20 LET D=1900 

225 REM R:AVERAGE RADIUS OF LATEX SPHERE 

230 LET R=5eE-07 

235 REM D9:PLATE SEPARATION 

2AG LET D9=eG2 

245 REM Pl:VALUE OF 'PI°* 

250 LET P1=3.14159 

255 REM MtAVERAGE MASS OF LATEX SPHERE 

260 LET M=4*P1*R*R*R*D/3_ 

265 REM Ki»K2,s AND K3 ARE USEFUL CONSTANTS 

270 LET K1=Q/D9 

275 LET K2=G*D9O*xM 

288 LET K3=6*P1+*H 

285 LET NO=@ = 
290 PRINT "INSTRUCTIONS (1=YES» @=NO)"3 = 
295 INPUT I : 
300 IF I=@ THEN 360 | 
365 IF I <> 1 THEN 299 | " 
319 PRINT | = : 
315 PRINT “INSTRUCTIONS -- AFTER EACH QUESTION MARKs"'3 

320 PRINT " (V= ?)5 YOU Mays" 

325 PRINT 

339 PRINT "TYPE IN VOLTAGE BETWEEN -199@ AND 1090 CIN ORDER TO MAKE" 
335 PRINT "THE VELOCITY PRINTED OUT AS CLOSE TO ZERO AS POSSIRBLE)»" 
349 PRINT “REQUEST CALCULATION OF CHARGE FOR STOPPED DROP"'s 

345 PRINT " (TYPE IN 2000).5"3 : 
356 PRINT "REQUEST NEW BATCH OF DROPS CTYPE IN 3009).5" 2 
355 PRINT "OR END THE PROGRAM (TYPE IN 4000)."" | 














¢ 360 PRINT 
365 PRINT 
379 PRINT " 'y"NQ ELECTRIC FIELD" 
375 PRINT 
380 PRINT “DROP: "sN3+1sNO+2sN0+3eNG+4 
385 PRINT e¢ FU FE ce ae we FE EF ce me ae OF, TPs we ae OF, OF ee OF 
390 IF N@>@ THEN 415 
395 PRINT " VELOCITY" 
40% PRINT "(METERS/SEC)" 
> 405 PRINT "¢€ X 101-6 9" 
419 REM GENERATE RANDOM CHARGE, RADIUS» AND MASS 
415 FOR I=1 TO 4 
Q 429 LET NCII=INTC3+3*COS(6-283*RNDC1))*SQRC-2*LOGCRND(1)))+605) 
425 IF ABSC(NCI)-4)>4 THEN 420 
430 LET RCI)D=R*Ce994+eG2*RND(I)) 
435 LET MCI)=4*P1*RCT)*RCID*RCI*D/3 
449 NEXT I 
445 PRINT " *, 
450 FOR I=1 TO 4 
455 PRINT e1*INTC1e E+07*% (KI#VENCI)=-MCI)*G)/(CK3#RO1))405)>5 
469 NEXT I 
465 PRINT" " 
470 PRINT "Vs"syst 3 
475 INPUT I 
489 IF ABSCI)>100@ THEN 495 
485 LET vsI 
499 GOTO 445 
495 PRINT 
569 IF I <> 2900 THEN 565 
505 PRINT 
518 PRINT “CALCULATION FOR WHICH DROP"; 
515 INPUT I 
520 FOR J=1 TO 4 
525 IF N@+tJ-I=9 THEN 540 
530 NEXT J 
535 GOTO 595 
S54 PRINT | 
545 PRINT "CHARGE ON DROP"SI3" 1S'3 -G1XINTC1eE+21*K2/VF05)3 
* 550 PRINT " X 19t-19 COULOMRBS." 
+ 555 PRINT | 
566 GOTO 479 
| | 565 IF I <> 3000 THEN 580 
a 570 LET NO=N0+4 
575 GUTO 375 
589 IF I <> AGOGO THEN 470 
585 END 

















ES a Se eee ee 
APPENDIX D. SAMPLE RUNS 
NO ELECTxIC FIELD 
DROP s 1 2 3 e: 
VELOUITY 
(METEKS/SEC) 
CZA l0t-6 2d 
=30«5 “2YeT7 = 3U ~3Ve 7 
V= O ? 1L0U 
 =256t -15e4 -2U e5 ~30e7 
V= 100 7° SOU 
et Aled 174 -~30e7 
V= 500 ? 6UU 
-2e3 56 2609 -30e7 
V= 600 t 6ea5d 
~lel 5925 2Ye3 -30e7 
V= 625 £650 
el 63el : 31-6 -3Ve7 
V= 650 ? 649 
O 63 3i<5 -30¢e7 
V= 649 ? 2VUVvvV 
CALCULATION FOr WHICH DROP? 1 
CHARGE ON DkOF 1 IS 1-56 xXx 10t-19 COULOMBS. 
V= 649 ? 2UVU 
-2lel ~l.el -11 -30¢7 
V= 2Uv 240 
-2U 06 e3 -10e1 -30-7 
V= 210 ? 2VU5 
-2U0 e7 U =lUsea = ~3UVe 7 
V= CUS 2? 20U0V0 
CALCULATION FOR WHICH DOr? @ 
CHARGE ON DROP 2 IS 4-93 aA 10t-19 COULOMBS-. 
8 : 





FRE Ee CE EES RR Ee IEEE OE FIERA NS RRR 5S Ar aT PCy ETE Pat NE 9 tee 





2 & Ve 208 ? 400 : 
2 | “1107 2704 79 “ses 
V= 400 ? 350 
-i4 20 e3 3e2 -30 e7 


V= 350 7? 325 


-15e2 l67 os - 3067 
‘t& ve 325 7 320 
ao -15-4 16 4 -30 e7 
e e320 2-315 
: -~15e7 1563 -~e] -30¢7 
Ve: 313: S20 
~15¢6 15-4 UO -30 e7 


V= 316 ? 2000 


CALCULATION FORK WHICH DROP? 3 : 
CHARGE ON DROP 3 I5 3-25 A 1Ut-19 COULOMBS-. 


V= 316 ? 3000 





DROP ro) 6 7 G 
74 43el =Loso 59 ed 


V= 316 ? 100 





Coed = Jon -~25¢6 “15 
V= 100 ? 90 

-06 “Sua -26e¢l -404 
= 9U 2? 92 

el “Ged -26 -~30d 
Vea yen FF 9] 

-e2 ; “9e5 -26 ~Ael 
Ye 94.7 91.7 ; 

U -“9e3 -26 ~“3eY 


V= GQle7T ? evUvU 


CALCULATION FO WHICH DROP? 5 





e CHARGE ON DHOP 5 IS 11219 « 101-19 COULOMBS.- 










































S43. 7 150 
3 1902 4-3 “233 lee7 
V= 150 ? 140 
1509 1.9 -2367 99 
V= 140 ? 135 
14.3 0 & -24 Ged 
Ve 135 2132 
13,5 i= =-2461 706 
V= 132 ? 131 
13 ~ 2 -2462 T03 
Ue 331 2 137<7 
13-2 6) -~24el 7:5 
V= 131-7 ? 2000 
CALCULATION FOR WHICH DROP? 6 
CHARGE ON DROP 6 1S 7679 & 10t-19 COULOMBS. : 
V= 131-7 ? 200 = 
3507 16 “20 09 2OeY = 
V= 200 ? 400 
101-7 620¢7 -11.e-5 bB3eS 
V= 400 ? 6U0 
167¢7 109 e5 -2el 140-7 
Ve 600 ? 625 
176 11655 =. 147% 
V= 625 ? 630 
177+6 116¢5 -e if 149 e2 
V= 630 ? 650 
16402 121.61 3 15409 
V= 650 ? 640 
140 9 11668 so 152 
V= 640 ? 642 : é 
1616 11923 -el 1526066 oe 
Ve 640-2 fas... 
. 1619 11965 =ei 15269 
V= 643 ? 644 
162-2 119-7 et 15342 
V= 644 ? 2000 
CALCULATION FOR WHICH DHOP? 7)” 
| CHARGE ON DROP 7 IS 1-59 X 10t-19 COULOMBS. 


















‘ 
EB 


1.3 
404 “602 “2504 -el 


4e7 => —-25e3 ea 


V= 106 ? 150523 
4e5 -~602 “25-4 U 


V= 105e3 ? 2000 


CALCULATION FOR WHICH DROP? & 
CHAHGE ON DROP & IS 9e75 A 1Ut-19 COULOMBS. 


V= 105e3 ? 4000 
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APPENDIX E. Falling Spheres 


Layer 1: 


ve 





“XK 


If a tiny sphere falls in vacuum under the action of gravity alone, the 
velocity at time t after it starts to fall is 


v=-gt. 


A graph of velocity vs. time looks like: 


V 


e 





“Vv 


If the same tiny sphere falls in air under the action of gravity, the 
velocity at time t is: | 


V; is a constant (the terminal velocity). A graph of velocity vs. time 
looks like: 





12 

















The terminal velocity is not reached as soon as the sphere starts to 
fall, but for the tiny latex spheres in this experiment it takes only 
about ten microseconds to reach the terminal velocity. The terminal 
velocity of a tiny sphere can be calculated: 


sof aes -mg/K 
where m is the mass of the sphere, g is the acceleration due to gravity, 
and K has to do with the drag caused by the air. If the radius of the 
sphere and the viscosity of air are known, K can be calculated. For 
the spheres used in this experiment and air at room temperature and 
pressure, 


gw ee newton-second/meter 


Since the mass of these spheres is 5.2 x io kg., the terminal 
velocity is about 2 millimeters per minute. 


When the electric field is turned on, the terminal velocity of a 
charged sphere becomes: 


v, * (ong Fel) ce 
where c is the charge on the sphere and E is the value of the electric 
field. E can be calculated from the voltage across the plates and the 


Spacing between them by, 


E=V/D. 


The graph of velocity vs. time looks like: 


V 


When the upward force due to the electric field is jist equal to the 
downward force due to gravity, the velocity is zero and the sphere 
stops. 
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Layer 2: 





‘© 





—X 


The equation for the motion of the charged sphere is 
F = ma = mv 


F is the total force on the sphere. It is made up of: 


F 


g force of gravity = -mg 


F, = bouyant force (since the sphere "floats" 
+0 Sit) * ts V8> 


where m,;, is the mass of air displaced by the sphere; 


& 
| 


Fy = drag force on the sphere due to its motion 
through the air = -Kv, 


where K can be calculated from the radius of the sphere and the viscosity 
of air; 


F, = electric force on the charge due to the electric 
field = cE. 


The bouyant force is less than one percent of the force of gravity and will 
be ignored here. It must be taken into account in precise determinations 


of the electronic charge. 


The equation can be written: 





mv = F = mg + cE - Kv 


This differential equation involving both the veanesty and the time rate 
of change of velocity has the solution: 


“t/to) 


v= v,~(1 -e 




















This is the equation that is plotted in the second graph on the first 
page. The velocity starts at zero and approaches v. at times long 
compared with t,. (When t = 3ty, v is already 95% of V+) Since 
for this form of v : 


: ~. -t/t 
Vv = : e / O 


v and v can be put back into the equation of motion, and the equation will 
be true if 


< 
I 


(-mg + cE)/K 


and 


t m/K . 


O 
The calculation of K is a problem in fluid dynamics. In 1843, Stokes 
found K, the drag force on a sphere falling slowly in a viscous fluid 
to be: 


K = 677hR 


where R is the radius of the sphere and h is the viscosity of air. The gen- 
eral problem of a body moving through a viscous fluid is very complicated 
and the sphere is one of the few shapes for which K can be calculated 

at all. Luckily, both oil drops and latex particles are very good spheres. 
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